Disorder-induced cubic phase in Fe2-based Heusler alloys. 
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Based on first-principles electronic structure calculations, we analyze the chemical and magnetic 
mechanisms stabilizing the cubic phase in Fe2-based Heusler materials, which were previously pre- 
dicted to be tetragonal when being chemically ordered. In agreement with recent experimental data, 
we found that these compounds crystallize within the so-called "inverted" cubic Heusler structure 
perturbed by a certain portion of the intrinsic chemical disorder. Understanding these mechanisms 
is a necessary step to guide towards the successful future synthesis of the stable Fe2-based tetrag- 
onal phases, which are very promising candidates for the fabrication of rare-earth-free permanent 
magnets. 
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One of the oldest problems within the field of ma- 
terials science is the search for the inexpensive hard 
magnets, i.e., for materials retaining their magnetisa- 
tion after being once magnetized. Their role in the 
daily life can be hardly overestimated: hard magnets are 
widely used in automotive applications, telecommunica- 
tions, data processing, consumer electronics, instrumen- 
tation, aerospace and bio-surgical applications. In par- 
ticular, they play a unique role in renewable energy tech- 
nologies based on electric generators (e. g., rotors in wind- 
turbines, small hydroelectric systems etc.). However, ma- 
terials exhibiting outstanding hard-magnetic properties 
together with high magnetization and high Curie tem- 
perature are rather expensive as being based on combina- 
tions involving rare-earth elements (e.g., Sm-Co, Nd-Fe- 
B) [l|,[2|. Thus, the development of new inexpensive com- 
pounds with hard- magnetic properties (i.e., rare-earth- 
free hard magnets) which can be industrially mass pro- 
duced is important and highly relevant (for the review 
see e.g., Ref. p|). The recent explosion of attention for 
the tetragonally-distorted magnetic Heusler systems orig- 
inates at a large extent from this prospective as well [3|- 

Indeed, apart of being promising candidates for tun- 
neling magneto-resistance and spin-torque-transfer ap- 
plications ll, U , this family may also provide materials 
combining the tetragonal distortion with a large mag- 
netic moment and high Curie temperature, suitable as 
hard magnets. The group of Fe2YZ-based Heusler com- 
pounds (with Y and Z being the transition and the main- 
group element, respectively), theoretically predicted to 
be tetragonal with a large magnetization (4-5 /ie/f-u., 
f.u.=formula unit) would be one of such promising ma- 
terials sources [5]. In contrary, the subsequent synthesis, 
XRD and Mossbauer characterization have shown that 
all these compounds crystallize in the cubic phase [gJ . To 
understand which ingredients can lead to their tetrago- 
nal distortion obviously implies an important preliminary 
step - a detailed understanding of the mechanisms sta- 
bilizing their cubic phase. This is the main point of the 
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FIG. 1. (color online) The cubic unit cell of the point- 
symmetry group No. 216 together with its schematic graph- 
ical diagram and the corresponding written notation. The 
high-symmetric Wyckoff positions 4a, 46, 4c and 4d are dis- 
tinguished by different colors. 



present study. 

Before proceeding to the results, first we would like 
to introduce the notations extensively used throughout 
in the text (see Fig. [1]). In the most general case any 
cubic Heusler system corresponds at least to the point- 
symmetry group No. 216. In order to distinguish between 
different chemical configurations, we will use the special 
written notation according to the occupations of the four 
high-symmetric Wyckoff positions: first we will write 
down the occupants of the 4a and 46 sites followed by 
the slash sign, then - of the 4c and 4(i, i. e., "4a46/4c4(i". 
Thus, e. g. the so-called "regular" and "inverted" vari- 
ants (terms introduced in Ref. [5]) of Fe2CuGa can be 
written as CuGa/FeFe and FeCa/CuFe, respectively. In 
case if certain Wyckoff position is occupied by several 
atomic sorts randomly, e.g. by A and B with probabil- 
ities X and 1 — X, it is noted using the square brackets: 
[Aa;Bi_a;]. In case of the tetragonal distortion (c ^ a) the 
symmetry reduces at least to the point group No. 119. 
The sequence for the written notation in this case does 
not change, it implies only the usage of the different 
Wyckoff positions: 2a2b/2c2d. 

To clarify the discrepancy between the theoretical pre- 
dictions from Ref. [5|] and the experimental data [6], in 
the following we will study the relative stability of the 
cubic and tetragonal phases of the Fe2-based Heusler 
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systems by optimizing both structural, magnetic and 
chemical degrees of freedom based on ab-initio density- 
functional calculations. As a suitable numerical tool 
which accounts for these factors simultaneously, we use 
the fully-relativistic Green's function formalism imple- 
mented within the SPR-KKR (Spin-Polarized fully Rel- 
ativistic Korringa-Kohn-Rostoker) method Q. The ran- 
dom occupation is described in terms of the CPA (Coher- 
ent Potential Approximation) [8|. Despite its mean-field 
nature (the effective averaging of the short-range order 
effects) CPA remains the most practical technique which 
includes the essential features of randomness. In order 
to ensure that the CPA result is not an artifact of the 
single-site approximation, we performed additional su- 
percell calculations. It was also found, that the usage 
of the full potential (i.e., the non-spherical potential) is 
much more essential for the adequate description than 
a particular choice of the exchange-correlation potential. 
For this reason, the presented calculations corresponds 
to the fully-relativistic and full-potential results employ- 
ing the local density approximation for the exchange- 
correlation functional [9|. The calculations for different 
c/a ratios are performed for the fixed volume, which was 
taken from the available experimental data |6i]. 

In order to explain the mechanism which keeps the 
Fe2-based systems cubic, throughout the discussion we 
will focus on the Fe2CuGa system, because as we found, 
all basic conclusions valid for this system can be trans- 
ferred without restrictions onto other compounds in this 
series (i.e., Fe2CuAl, Fe2NiGa, Fe2NiGe and Fe2CoGe) 
synthesized experimentally [6[. 

The main outcome of the present study is summarized 
in Fig. [21 which represents the dependency of the total 
energy on c/a ratio for various alloy configurations. The 
ordered "regular" Heusler structure in cubic phase turns 
out to be unstable (indicated by the corresponding en- 
ergy curve maximum at c/a = 1, black line in Fig.[2](a)), 
whereas at about c/a = 1.54 the system falls into the 
relatively deep energy minimum. For the fixed chem- 
ical order (i.e., CuGa/FeFe) the tetragonal distortion 
is the only mechanism which can relax the instability 
of the cubic phase, since the magnetic degrees of free- 
dom are already in use (for more detailed description of 
tetragonal distortion mechanisms, see e. g. Si). As 
it follows, the gradual transition towards the "inverted" 
Heusler structure, realized by random chemical Cu-Fe 
inter-layer exchange starts gradually to develop an en- 
ergy minimum for the cubic phase. Although the config- 
urations with intermediate Cu-Fe site occupations (e.g., 
[Cuo.5Feo.5]Ga/[Cuo.5Feo.5]Fe) exhibit an energy minima 
for both tetragonal and cubic phases, the limiting ordered 
system {x = 1, i.e. the fully "inverted" FeGa/CuFe) is 
stable only within the cubic phase (Fig. [2] (a) or (b), 
red curve). Despite the large energy difference (about 
— 170 meV/f.u.) between the cubic "regular" and "in- 
verted" phases, the deepest absolute energy minimum 
is found for the tetragonally-distorted "regular" config- 
uration (see Fig. [2](a)), which is by about 20 meV/f.u. 
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FIG. 2. (color online) Total energy of Fe2CuGa Heusler al- 
loy calculated as a function of c/a ratio for various distri- 
butions of Fe and Cu (indicated by the box-like diagram; 
green, yellow and white colored areas correspond to the Fe, Cu 
and Ga occupations, respectively), (a) Thick black and red 
curves correspond to CuGa/FeFe ("regular") and FeGa/CuFe 
("inverted") configurations, respectively. Thinner violet line 
shows the intermediate [Cuo.5Feo.5]Ga/[Cuo.5Feo.5]Fe case be- 
tween "regular" and "inverted". The absolute energy min- 
imum of the "regular" (CuGa/FeFe) tetragonal phase is 
taken as a reference. (b) Solid red curve represents the 
same "inverted" configuration as on the left panel, whereas 
the solid blue line shows the most stable configuration: 
FeGa/[Cuo.5Feo.5] [Cuo.sFeo.s], obtained from the ordered "in- 
verted" by mixing Fe and Cu randomly in-plane. The dashed 
red and blue lines correspond to the distributions derived from 
the previous two configurations via additional in-plane ran- 
dom spread of Ga and Fe: [Feo.sGao.s] [Feo.sGao.sJ/CuFe (red 
dashed) and [Feo.sGao.s] [Feo.sGao.sJ/iCuo.sFeo.s] [Cuo.sFeo.s] 
(blue dashed). The absolute energy minimum of the "in- 
verted" (FeGa/CuFe) phase is taken as a reference. 



more stable compared to the "inverted" cubic configura- 
tion. Thus, for the ordered systems our results agrees 
with the former calculations [5|. This means, that the 
mechanisms stabilizing the cubic phase involve degrees 
of freedom which where neglected so far, e. g. the chem- 
ical disorder [6|]. 

It is important to note, that the huge energy decrease 
(about —170 meV/f.u.) gained by going from the "reg- 
ular" cubic to the "inverted" cubic configuration (the 
largest energy scale in the diagram on Fig. [2j) is most 
likely of the magnetic origin. The latter is due to the op- 
timization of the magnetic exchange coupling within the 
Fe sublattice, since the nearest magnetic neighbors (i.e., 
the Fe atoms from the adjacent layers within "inverted" 
cubic FeGa/CuFe) are sitting closer to one another com- 
pared to the "regular" cubic CuGa/FeFe setup, in which 
they are in-plane. For this reason, by searching for the 
more stable configurations, we start from the "inverse" 
cubic system and perturb it by in-plane chemical disor- 
der (i.e., by conserving the total amount of Cu and Ga 
within adjacent layers). Hence, there are two important 
in-plane disorder scenarios: random in-plane mixtures of 
Fe-Ga, and that of Fe-Cu. 



3 



We found, that the random in-plane spread of Ga and 
Fe (case [Feo.5Gao.5][Feo.5Gao.5]/CuFe, red dashed hne 
in Fig. [21(b)) leads to the increase of the total energy 
(compared to the "inverted" configuration, FeGa/CuFe) 
by about 150 meV. In contrast, the random in-plane 
spread of Cu and Fe (case FeGa/[Cuo.5Feo.5] [Cuo.sFeo.s], 
solid blue line in Fig. [2](b)) leads to an energy gain of 
about —40 meV (again, compared to the "inverted" case, 
FeGa/CuFe). The key observation is, that this —40 meV 
energy gain is enough to stabilize the cubic structure 
(in FeGa/[Cuo.5Feo.5] [Cuo.sFeo.s configuration), which fi- 
nally becomes more stable than the tetragonal "regular" 
ordered CuGa/FeFe by about 40 - 20 = 20 meV/f.u., in 
agreement with experiment. 

Our results show, that these two effects (Fe-Ga and Fe- 
Cu in-plane random mixtures) are rather independent on 
one another: i. e. disregarding the particular arrangement 
of atoms within the adjacent layer, the energy changed by 
150, 40 or 150 ± 40 meV while going from one distribu- 
tion to another within all four cases. The large increase 
in energy by 150 meV in the first case is mainly due to 
the distinct nature of Fe and Ga. So, within the fixed 
square lattice it is unfavorable to form separate clusters 
of Fe and Ga, since each atomic sort would prefer to cre- 
ate its own lattice within a cluster which will be rather 
different from another. For this reason, any perturbation 
of the perfect chemical order in Fe-Ga layers will increase 
the total energy. This issue, however, is not critical for 
the second case: the separation of Fe and Cu within the 
given lattice does not cost so much energy, since both 
atom types are much more similar. To ensure that the 
—40 meV energy gain in this case is not just an artifact of 
the single-site nature of the CPA, we have performed su- 
percell calculations by systematically increasing the num- 
ber of Fe-Cu in-plane swaps, mimiquing Fe-Cu disorder. 
This has shown that by increasing the degree of Fe-Cu 
separation the total energy is indeed reduced by around 
-40 meV. 

The subsequent calculations of the magnetic exchange 
coupling constants Jij (Fig. [3]) of the classical Heisenberg 
model {H = — Xli>j Jij^i^j-> where e^^j are the unity vec- 
tors along the magnetization directions on local sites i 
and j) revealed the magnetic origin of both stabilization 
mechanisms responsible for the atomic rearrangement 
from the "regular" into the "inverted" phase and for the 
chemical disorder within the Fe-Cu layers. Namely, the 
strong Fe-Fe inter-layer coupling ( Jinter-iayer ^ 25 meV 
between the adjacent Fe-Ga and Fe-Cu layers) keeps the 
whole system ferromagnetic. This is in agreement with 
the high Curie temperature (798 K) measured in Ref. 
Although the in-plane couplings appear to be an or- 
der of magnitude weaker, still, as we mentioned, their 
optimization plays a crucial role in the stabilization of 
the cubic phase. In the ordered "inverted" configura- 
tion (FeGa/FeCu) the nearest in-plane Fe atoms tend 
to couple antiferromagnetically (Jin-piane = — 1.4 meV). 
This interaction works against the overwhelming fer- 
romagnetic order already set by the strong inter-layer 




FIG. 3. (color online) Comparison of the magnetic exchange 
coupling in the "inverted" FeGa/FeCu (a) and in the most 
stable FeGa/peo.sCuo.s] [Feo.sCuo.s] (b) configurations. The 
atoms are arranged within Fe-Ga and Fe-Cu layers marked by 
light- and dark-blue horizontal planes, respectively. Fe, Cu 
and Ga atoms are shown as green, yellow and white spheres, 
respectively. Magnetic moments are shown by arrows. The 
bond thickness reflects the strength of the exchange inter- 
action. The inter-layer Fe-Fe interactions are dominating: 
^inter-layer ~ 25 mcV (thick grccu bouds) . The in-plane in- 
teractions are negligibly small, except those in Fe-Cu planes. 
Case (b) illustrates the typical distinction from the ordered 
"inverted" structure: the random in-plane swap of one Fe and 
one Cu atoms which brings two Fe atoms closer to one another 
within the Cu-Fe plane. This alters the nearest in-plane Fe-Fe 
exchange from antiferromagnetic (Jin-piane = — 1.4 meV, thin 
red bond) in case (a) to ferromagnetic (Jin-piane = 5.2 meV, 
thin green bond) in case (b). 



coupling. Thus, the magnetic energy can be further 
reduced by bringing the Fe atoms closer together as 
shown in Fig.[3](b), so that they couple ferromagnetically 
(^in-plane = 5.2 mcV). Put in practical terms, this effect 
favors the formation of random Fe clusters within the Fe- 
Cu planes, which -in contrast to the ordered case- can be 
more adequately described by the chemical disorder pic- 
ture, i.e. by implying the FeGa/[Feo.5Cuo.5] [Feo.sCuo.s] 
configuration. 

The considerations presented above can be futher sup- 
ported by comparing the electronic structures as shown 
on Fig. m The instability of the electronic subsystem 
is typically related to the strength of the DOS peaks in 
the vicinity of the Fermi energy. In case of the "regular" 
CuGa/FeFe cubic system, a huge instability peak at 
(total DOS(£^f) ^ 6.1 sts./eV) is produced by the van- 
Hove singularity in the minority-spin channel at W-point 
of the Brillouin zone (Fig. |4](a)). By replacing Cu from 
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FIG. 4. (color online) Comparison of the spin-resolved 
electronic band structures and related densities of states 
for the (a) "regular" cubic CuGa/FeFe, (b) "inverse" cubic 
FeGa/CuFe and (c) FeGa/peo.sCuo.s] [Feo.sCuo.s] configura- 
tions. The majority- and minority-spin states are distin- 
guished by red and blue, respectively. 



the 4a (or 4b) site will split the minority-spin states at 
W-point far away from the Fermi energy, noticeably re- 
ducing the DOS peaks (see Fig.|3](b) and (c)). This step 
is related to the largest ( — 170 meV/f.u.) energy gain at- 
tributed to the inter-layer magnetic exchange optimiza- 
tion discussed above. Now by going from the ordered 
"inverted" Heusler variant (FeGa/CuFe) to the most sta- 
ble case of FeGa/[Feo.5Cuo.5][Feo.5Cuo.5], the DOS at 
is further reduced (from 3.2 to 2.9 sts./eV). Indeed, 
as we have shown above, the Fe-Cu disordered configu- 



ration benefit from the in-plane magnetic optimization, 
and therefore it is by 40 meV/f.u. lower in energy com- 
pared to the "inverted" FeGa/CuFe system. We would 
like to point out, that very similar stabilization mech- 
anisms characterized by comparable energy scales can 
take place in other Fe2-based cubic Heusler compounds. 
For example, for Fe2CuAl and Fe2NiGe the inter-plane 
exchange energy optimization (i. e. by going from the 
"regular" to the "inverted" cubic phase) gains —367 and 
— 168 meV/f.u., whereas the in-plane optimization (due 
to Fe-Y in-plane disorder) contributes with —27 and 
—40 meV/f.u., respectively to the energy. Thus the sta- 
bilization mechanisms presented in this letter are rather 
general within the group of Fe2YZ materials. 

To conclude, we emphasize that the presented analysis 
explains the stability of the cubic phase in Fe2 YZ Heusler 
compounds, and provides a clear explantion for the dis- 
crepancy between experimental results and theory. The 
actual stabilizing mechanism appears to be the chemi- 
cal disorder, which optimizes the magnetic exchange cou- 
pling within Fe-Y layers of the initially ordered FeZ/FeY 
cubic phase. At the same time, the FeZ layers remain 
chemically ordered due to a large difference (i. e. atomic 
radius, valency, electronegativity etc.) between Fe and 
the main-group element Z. Thus, the most stable config- 
urations can be written as FeZ/[Feo.5Yo.5] [Feo.sYo.s]. In 
general, the important prerogative enabling the chemi- 
cal disorder is the "inverted" ordered structure: as we 
have seen, the effect of the rearrangement from YZ/FeFe 
into FeZ/YFe within the cubic phase is comparably ef- 
ficient to the tetragonal distortion in YZ/FeFe. On the 
other hand, the first scenario allows to further optimize 
the system by chemical disorder, whereas the second one 
does not. 
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